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Abstract. We present the first near-infrared (A^-band) homogeneous observations of a complete sub-sample of the 
3CR radio catalogue comprising all High Excitation Galaxies (HEGs) at z <0.3. After showing that the surface 
brightness decomposition technique to measure central point-like sources is affected by significant uncertainties for 
the objects in the studied sample, we present a new, more accurate method based on the R — K' color profile. Via 
this method we find a substantial nuclear A'-band excess in all but two HEGs - most likely directly associated to 
their nuclear emission - and we measure the corresponding 2.12 /im nuclear luminosities. Within the frame of the 
unification scheme for radio-loud active galactic nuclei, it appears that obscuration alone is not able to account 
for the different nuclear properties of the majority of the HEGs and Broad Line Radio Galaxies (BLRGs), and 
also scattering of the (optically) hidden nuclear light from a compact region must be invoked. More precisely, for 
~70% of the HEGs the observed point-like optical emission is dominated by the scattered component, while in the 
if '-band both scattered and direct light passing through the torus contribute to the observed nuclear luminosity. 
The estimated fraction of scattered light ranges from a few tenths to a few percent, while the torus extinction is 
between 15< Ay, torus <50 mag with only a few exceptions with lower obscuration. 
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1. Introduction 

Unification models for radio-loud active galactic nuclei 
(RL AGNs) aim at understanding the different observed 
properties of intrinsically identical so urces in terms of 
aniso tropics in the nuclear emission (see lUrrv fc Padovanil 
Il995l for a review). Anisotropy can arise from relativistic 
beaming and/or selective obscuration caused by material 
arranged in a non-spherical, torus-like geometry. 



RL AGNs can be broadly divided into two classes 
based on their extended radio morphology: low-power, 
edge-darkened FR I and high-power , edge-brightened 
FR II sources IIFanar off & Rilevl Il974lh For the power- 
ful FR II radio-galaxies (RGs), the unification of the 
broad-lined (BLRGs) and narrow-lined objects (NLRGs) 
invokes the existence of an obscuring 'torus' blocking the 
direct view of the nuclear and broad-line emission at 
large angles with respect to the jet axis. The detection 

(a still small num- 
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ber of) NLRGs (e. 


y. Antonucci 


1984: iGohen et alJll99<* 


lYounsr et alJl996a: 


Oele et alJl997l) provides the most di- 



rect support to this scenario. Conversely, in low-luminosity 
FR I RGs obscuration is actually neither 'requir ed' nor 
supported observationally l|Urrv fc Padovanilll995|) . 

A complication in this scheme arises due to the spec- 
tral inhomogeneity of the FR II RG population, compris- 
ing within the NLRGs class both high- excitation ( HEGs 



and low-excitation (LEG s) galaxies l|Laing et al.l Il99 
Ijackson fe Rawlingslll997|) . while all FR Is are character 
ized by low-excitation spectra. 

The study of the nuclear properties of RGs strongly 
benefited from the results of Hubble Space Telescope 
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(HST) optical and UV imaging of most 3CR RGs 
(jspinrad et all Il985|) with redshift z <0.3 (33 FR Is 
and 63 FR II; Chiaberge, Capetti & Celotti 1999, 
2002, CCC99, CCC02 hereafter; Chiaberge et al. 2002). 
Unresolved nuclear sources are detected in most of the 
FR Is, suggesting that their nuclei are essentially un- 
obscured (CCC99). The observed large range in the 
optical-UV spectral index can be accounted for just by 
moderate absorption (Ay ^1-2 mag), due to either ex- 
tended (kpc-scale) dus t lanes or (~100 pc) dusty disks 
l|Chiaberge et al.ll2002(h Furthermore, the nuclear optical 
luminosities strongly correlate with the 5 GHz radio core 
ones, arguing for a common synchrotron origin, and in 
turn correspond to extremely low Eddingt on ratios (e.g. 
Fabian fc Reeslll995l IZirbel fc Baumlll995l IWillott et all 
2000 HWills et al.l2004 . Thus FR I RGs might be lacking a 



l|Hill. Goodrich fc DePovlll996l) and thermal-IR (L'~ and 



substantial broad line region, torus and thermal disc emis- 
sion. Unresolved nuclear components are also detected in 
a significant number of FR lis (CCC02). With respect 
to FR Is they show a more complex behaviour which is 
however clearly related to their spectral classification: i) 
BLRGs have optical cores with luminosities in excess to 
those of FR I of similar radio power, possibly dominated 
by thermal emission; ii) no central source is seen in several 
HEGs, likely constituting the obscured narrow-lined coun- 
terparts of BLRGs; iii) in a few objects (e.g. 3CR 109 and 
3CR 234) the detected nuclear component is consistent 
with being the (obscured) quasar light scattered into our 
line-of-sight; iv) LEGs harbor faint optical cores, essen- 
tially indistinguishable from those seen in FR Is, indicat- 
ing that the FR I /II dichotomy is not univocally connected 
with the nuclear structure (LEGs constitute a substantial 
fraction of FR II sources). 

Nonetheless, the optical results can be strongly ham- 
pered by the effects of dust extinction. This problem is cer- 
tainly serious in RGs where (even leaving aside the pres- 
ence of nuclear tori) nuclear discs and dust lanes do not 
allow us to study the central regions in a significant frac - 
tion of the 3CR sources (~ 30 %) (Ide Koff et al.ll2000l) . 
Particularly instructive is the case of Centaurus A, the 
nearest RG. An unresolved nuclear compo nent is seen in 
HST images from the V to the K band I Marconi et, alJ 
|2000|) but its infrared flux is five orders of magnitude larger 
than the optical one. A similar effect is also commonly seen 
in Scyfert 2, where HST imaging have revealed the ubiq- 
uitous presence of nea r-infrared (NIR) nuc l ei in optically 
obscu red sources (e.g. iKulkarni et al.lfl998t lOuillen et alJ 
l200lh . Thus to gain a better understanding of the physics 
of RL AGNs it is crucial to obtain a more unbiased view of 
RGs nuclei, which can be provided by NIR observations. 

Relatively lit tle effort has been devoted to date to NIR 
imaging of RGs llLillv fc Longairlll982HTavlor et aill996t 
ISimnson. Ward fc Wallll200(i hereafter SWW00) and no 
systematic prog r am h as been carried out. Moreover in 
iLillv fc Longairl l|l982|) only 10"-12" aperture photome- 
try has been performed, leading to only about a dozen 
objects with NIR nuclear luminosity estimates, out of the 
complete 3CR RG sample with z <0.3. NIR spectroscopy 



M-band) imaging ( Simpson. Rawlings fc LacvU iggQ) for 
a total of about 15 RGs have revealed quasar nuclei ob- 
scured by Ay from - - 15 to >45 mag (SWW00). 

The NIR results can/ will be also compared to those 
in the mid-infrared (MIR) and in the far-infr ared (FIR). 
Indeed the recent study of 3CR sources by lHaas et "all 
(2004) provided evidence for the presence in high-power 
RGs of powerful AGN hidden by a dusty torus: steep- 
spectrum radio quasar, BLRGs and HEGs showed very 
similar IR spectral energy distributions, characterized by 
a FIR dust emission component, with relatively weaker 
MIR, NIR and optical emission in HEGs with respect to 
BLRGs. On the contrary, FR Is and LEGs were character- 
ized by a different spectral energy distributions, showing 
a dust emission bump in the MIR-FIR from dust cooler 
than in FR II, a bump peaking at 1 /xm due to the ellipti- 
cal hosts and a steep synchrotron spectrum in the cm-mm 
wavelength range. A comparison between 3CR FR Is and 
high-power FR lis in the MIR and FIR has been also 
recently performed by iMiiller et al.l l|2004h showing that 
FR Is are characterized by MIR and FIR dust luminosi- 
ties much lower than FR lis, consistent with the former 
hosting low-luminosity nuclei. 

To extend the NIR diagnostic to a large complete sam- 
ple we performed a AT'— band survey of all RG belonging 
to the 3CR sample at z < 0.3, aiming at: i) identifying 
optically obscured sources from the presence of a NIR ex- 
cess with respect to the extrapolation of the optical (HST) 
flux; ii) establishing the role of obscuration in FR I RGs; 
iii) exploring the nature of narrow lined FR lis, in par- 
ticular those with faint optical nuclei; iv) comparing the 
NIR properties of low and high luminosity radio sources; 
v) studying the properties of the nuclear regions of the 
sources affected by large scale dust obscuration. 

In this paper we present results obtained for the sub- 
sample of FR II RGs with an HEG optical spectrum. The 
outline is the following. In § El we describe the IR obser- 
vations and data reduction; a comparison between the IR 
and HST images aimed at finding obscured nuclei is pre- 
sented in § El In § E| and § El we focus on the results and 
discuss some of their consequences. Our conclusions are 
drawn in the final § El A zero cosmological constant uni- 
verse with Ho=75 km s" 1 Mpc~ x and qo=0.5 is adopted 
throughout the paper for comparison with previous works 
in this field. 

2. The sample, observations and data reduction 

2.1. The sample 

Among the 3CR RGs at z < 0.3, the HEG sub-sample 
co mprises 32 objects (CCC02 ). Following the definition 
bv lJackson fc Rawlingj l|l997|) . HEGs have been classified 
on the basis of their [OIII] equivalent width (larger than 
10 A), and/or the [Oil] /[OIII] ratio (smaller than 1). In 
Table ^ names, redshifts and radio data of the objects 
are reported, together with the optical [OIII] line and core 
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continuum (HST) luminosities, as taken from the litera- 
ture. Six objects of the sample have been excluded from 
the analysis: 3CR 63, 3CR 135 and 3CR 192 were not ob- 
served because of scheduling problems or clouds, 3CR 153 
had very bad seeing, 3CR 33 and 3CR 105 have no avail- 
able HST image (see S EHjl . 

2.2. TNG observations 

Observations have been obtained with the 3.6m Telescopio 
Nazionale Galileo (TNG), the Italian national facility lo- 
cated at La Palma Island (Spain) in two runs on Jul 8-12, 
2000 and on Feb 9-13, 2001. During the first run we used 
ARNICA (ARcetri Near Infrared CAmera), a 256x256 
pixel NICMOS 3 array with a pixel size of 0'.'35 and a 
field of view (FoV) of l'.5xl'.5, while for the second run 
NICS (Near Infrared Camera Spectrometer), a Rockwell 
1024x1024 HgCdTe Hawaii array, was already available. 
In its small field mode the pixel size is 0'.'13 for a total 
FoV of 2'. 2x2'. 2, while in its large field mode pixel size 
and FoV are doubled. 

Images were taken with the dithering mode with the 
source located in a pattern of 20 positions with a width of 
±30". All sources of this sub-sample are sufficiently small 
that it was possible to always keep the source within the 
images during the pattern (for both ARNICA and NICS). 
For NICS, the individual integration time was set to 1 min 
to achieve background limited images, i.e. a total integra- 
tion time of 20 min for each source (except 3CR 184.1 and 
3CR 234, for which it was 16 and 15 min, respectively). 
For ARNICA the individual integration time was set to 
48 sec, and the total ranged from 4.8 to 24 min depending 
on the source. We used the K' filter that has a central 
wavelength of 2.12 /im and a FWHM of 0.35 /zm. 

2.3. Data reduction 

The images taken with NICS were reduced with SNAP 
(Speedy Near- infrared data Automatic Pipeline), 
a pipeline for automatic reduction developed by 
F. Mannucci. SNAP performs a full reduction, 
with flat-fielding, sky subtraction, computation of 
the offsets, correction for geometrical distortion, 
object masking and correction for cross-talk be- 
tween the quadrants (for more information, see 
http : //www. arcetri . astro . it/~f ilippo/ snap/|. 
For the ARNICA data, a flat field was produced by 
median filtering of the frames, after scaling each frame 
to have the same median pixel value. Each frame was 
then divided by a normalized flat-field and registered 
using as reference the target itself and, when possible, 
stars present in the FoV. The aligned images were then 
averaged to produce the final reduced image. 

For both NICS and ARNICA, absolute flux calibration 
was obtained observing 6-11 standard fields (depending on 
the night) using a 5 position dithering pattern. From each 
frame of each standard star we derived a measurement of 



the zero-point, for a total of 32-61 independent values per 
night. The dispersion of the zero-point measurements pro- 
vides an estimate of the calibration accuracy in the range 
4-6% for NICS and 3-6% for ARNICA, depending on the 
night. NICS last night was not photometric, with thin cir- 
rus, and no standard star has been observed. The absolute 
flux calibration was obtained using the stars present in the 
fields in common with previous nights (these fields were 
re-observed since the seeing in last night was much bet- 
ter). The flux calibration proved to be stable throughout 
the night, consistent with previous night calibrations and 
with an accuracy of ~ 6%. 

The absolute flux calibration has been compared with 
data published in the literature, available for 3CR 79, 
3CR 98, 3CR 1 71, 3CR 223, 3CR 23 4 (SWW00), 3CR 456 
and 3CR 460 (|de Vries et al.lll998l) . After correcting for 
the different filter used \K' = K + (0.22 ± 0.03)(tf - 
iO: lNeumann. Meisenheimer fe Roserll9 97], our aperture 
photometry is perfectly consistent with theirs (the mean 
difference is 0.02±0.03 mag, less than our absolute flux 
calibration error). The images of all sources - together 
with their corresponding HST ones, see S 13.21 - are shown 
in Appendix ITJ1 

3. Analysis 

3.1. The surface brightness decomposition technique 

The 'standard' technique used to detect and measure point 
sources at galaxy centers requires a one (or two) dimen- 
sional synthetic modeling of the galaxy profile to which 
a central point source is added. Setting the seeing from 
observations, free parameters of this procedure are those 
describing the galaxy profile (e.g. central surface bright- 
ness and effective radius) and the intensity of the nuclear 
source. This technique works well with HST images, not 
affected by seeing effects, and high quality ground-based 
photometry of local sources. Indeed only the HST spatial 
resolu tion and sensitivity to high surface brightness fea- 
tures l)Hutchings fc Morris! Tl 9 951: iMcLeod fc Riekelll995[) 
have allowed the detect ion of faint point-like sou rces at the 
centre of typical RGs l|Capetti fc Celottilll999t CCC99), 
totally swamped by the host stellar emission at ground- 
base resolution. 

We have applied this procedure to measure the IR 
nuclear luminosity Lk' of our sources, by using the two- 
dimensional fitting algorithm GALFIT l|Peng et alJl2002h 
designed to extract structural components from galaxy im- 
ages by simultaneously fit an arbitrary number of them. 
In Appendix [S] the results from such analysis are pre- 
sented and discussed for representative sources. In sum- 
mary we find that: 1) for the high z sources and those 
with faint nuclei (such as 3CR 405), the measured Lk' 
can differ by more than a factor of 20 depending on the 
initial parameters; 2) in many cases the fits are equally 
good with or without a central point source [even in the 
case of 3CR 405, in which this is clearly detected both in 
the HST image and with the R — K' color profile analysis 
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Table 1. The sample of HEGs 



Source Name 


Redshift 


Images" 


log Ll78 


log L r 


log LroTTTl 


log Lo 




2 




erg s _1 Hz -1 


erg s _1 Hz -1 


erg s _1 


erg s^ 1 Hz -1 


3C 18 


0.188 


H,A 


34.16 


31.84 


41.88 


28.62 


3C 33 


0.059 


A 


33.56 


30.27 


42.01 


not obs. 


3C 63 


0.175 


H 


34.08 


30.97 




28.29 


3C 79 


0.256 


H,N 


34.54 


31.19 


42.43 


28.33 


3C 98 


0.030 


H,N 


32.79 


29.25 


40.87 


<25.67 


3C 105 


0.089 


N 


33.45 


30.36 


40.78 


not obs. 


3C 135 


0.127 


H 


33.78 


30.18 




27.36 


3C 153 


0.277 


H,N 


34.44 


29.73 




<26.66 


3C 171 


0.238 


H,N 


34.42 


30.36 


42.50 


26.50 


3C 184.1 


0.118 


H,N 


33.61 


30.25 


42.19 


28.17 


3C 192 


0.060 


H 


33.22 


29.73 


41.60 


<26.99 


3C 197.1 


0.131 


H,N 


33.45 


30.30 


41. 26 6 


28.10 


3C 198 


0.082 


H,N 


33.38 




41.04 


28.03 


3C 223 


0.137 


H,N 


33.74 


30.57 


42.17 


<27.27 


3C 223.1 


0.108 


H,N 


33.31 


30.24 


41.65 


<27.18 


3C 234 


0.185 


H,N 


34.27 


31.88 


43.17 


29.10 


3C 284 


0.239 


H,N 


33.87 


30.25 


42.19 


complex 


3C 285 


0.079 


H,N 


32.87 


29.93 


40.73 


25.65 


3C 300 


0.270 


H,N 


34.45 


31.07 


42.16 


27.79 


3C 321 


0.096 


H,A 


33.32 


30.78 


42.15 


complex 


3C 327 


0.104 


H,A 


33.97 


30.88 


42.05 


<26.55 


3C 349 


0.205 


H,A 


34.07 


31.18 


41.42 


28.16 


3C 357 


0.167 


H,A 


33.67 


30.48 


41.76 c 


<26.90 


3C 379.1 


0.256 


H,A 


34.03 






<27.21 


3C 381 


0.161 


H,A 


33.91 


30.48 


42.37 


<27.52 


3C 402 


0.025 


H,A 


32.04 


29.73 




26.59 


3C 403 


0.059 


H,A 


33.28 


29.87 


41.55 


26.65 


3C 405 


0.056 


H,A 


35.73 


32.20 




complex? 


3C 436 


0.215 


H,A 


34.18 


31.21 


41.52 


complex 


3C 452 


0.081 


H,A 


33.88 


31.24 


41.35 d 


<26.89 


3C 456 


0.233 


H,A 


34.21 


31.39 


42.59 


28.50 


3C 460 


0.268 


H,A 


34.10 


31.39 


41.67 


27.60 


"Images coding: 


H=HST, A: 


= Arnica, N 


=NICS; 









L178 is the total radio luminosity at 178 MHz; L r is the nucl ear radio luminosity at 5 GHz (references in CCC02); 
L[om] is the luminosity of the [OIII] 5007 A emission line from I.Tackson fc Rawlingsl 1^997) and C. Willott's web page, 
http://www-astro.physics.ox.ac.uk /~ciw/3crr/3c rr.html; except for b 3CR 197.1 (from the SDSS spectrum), C 3CR 357 
iCrawford. Fabian fc Johnstonelll98 8) and d 3CR 452 llRawlings et alJll98^ : Lo is the luminosity of the optical core as mea- 
sured in CCC02. 



(see Appendix[BJ)]; 3) consistent and robust results are ob- 
tained only for those objects in which the nuclear compo- 
nent clearly dominates over the host galaxy (e.g. 3CR 456 
and 3CR 234). As this method does not provide sound 
results for objects with faint nuclei or at medium/high 
z, which represent the majority of our sources, we devel- 
oped and adopted a new procedure, presented in the next 
section. 

3.2. Comparison with HST observations: the color 
profile technique 

The alternative approach we developed takes advantage 
of the fact that HST images are available for the targets 
and thus the galaxy profile and (when detected) the in- 
tensity of the optical nuclear source are known a priori. 
Therefore, instead of modeling directly the IR images we 



compared them to the HST ones (F702W R band filter). 
The basic idea is that an obscured nucleus would reveal it- 
self through an increase of the IR nuclear flux with respect 
to what is seen in the optical images. 



In order to perform the comparison it is necessary to 
produce a synthetic HST image matching both the seeing 
and the pixel size of each IR image. Therefore we fitted a 
Gaussian to all the stellar objects in the HST field, mea- 
sured their FWHM and adopted as FWHM the median 
of all measurements with an uncertainty estimated from 
their dispersion. The values obtained range between 0'.'6 
and l'/7 with a typical uncertainty of ±7% (the error on 
the seeing is between 0"02 and 0"2). The HST image has 
been then convolved with a Gaussian with the appropri- 
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ate FWHM and interpolated to match the pixel size of 
NICS/ ARNICA 1 . 

The color-profile method is detailed in Appendix iBl for 
two show-case galaxies, namely 3CR 403 and 3CR 300. 
Briefly, photometry is performed on annuli of increasing 
radius out to 3" on the 'matched' HST and on the re- 
duced TNG images. An average color (and its associated 
error) is derived at each annulus. The outer portion of each 
color profile sets the color of the host galaxy. The nuclear 
monochromatic if'-band luminosity Lk' is estimated as 
the sum of Lk',xsj i-e. the measured nuclear IR excess, and 
Lk\ scaled; namely the HST optical nuclear component Lo 
properly scaled with the reference host galaxy color. 

Also for all RGs (but two) with an optical nuclear up- 
per limit, a significant IR excess is present in the inner 
part of the color profile. Clearly in such cases Lk> is less 
accurate: Lk', scaled is calculated considering Lo as a de- 
tection, and its error takes into account the possibility of 
Lq equal to zero. For 3CR 327 and 3CR 379.1 instead 
only an upper limit to the nuclear IR excess could be es- 
timated and therefore Lk' is also an upper limit. Four ob- 
jects (namely 3CR 284, 3CR 321, 3CR 405 and 3CR 436) 
have optical complex nuclei. In these cases, although no 
measurement is available for Lk', scaled - implying that the 
resulting Lk* should be strictly considered as a lower limit 
- the 'true' value would be very similar to the one esti- 
mated if Lk 1 , scaled <SC £k',xs- 

3.2.1. Host galaxy properties 

In Fig. n we present the color profile derived for each 
source. For the majority of them the host galaxy color 
is consistent with being constant; only a very few objects 
show significant gradients or complex behaviour in the 
R — K' profile (see Appendix [C] for details on individual 
objects). For 19/26 the outer part of the profile has been 
fitted with a linear model R — K'\ modci — A + Br in order 
to account for a possible gradient: 5/19 sources (3CR 98, 
3CR 285, 3CR 349, 3CR 357 and 3CR 405) show signifi- 
cant (at Jj3er) gradient, while other 7/19 objects (3CR 79, 
3CR 223, 3CR 300, 3CR 379.1, 3CR 403, 3CR 436 and 
3CR 452) have only marginally significant (>lcr) color 
gradient. The mean R — K' color of the host galaxy is 
2.97±0.29, consistent with the value of R- Jf=2.75±0.45 
of a sample of RGs and RL QSOs hosts ijDunlop et alJ 
120031) and with th e expected color R — K=2.8A for an el- 
liptical at z=0.16 l(Poggiantilll997h (the mean redshift of 
the sample is 0.15). The mean value of the R — K 1 gradient 
[defined as A(R- K')j A(logr)] is -0.22±0.02 mag dex" 1 , 
consistent with the values of: i) -0.14 mag dex -1 (with a 



1 We chose to use a Gaussian rather than a real PSF to con- 
volve the HST images because the Gaussian model reproduced 
well the observed PSF. A small dependence of the PSF shape 
on the position in the FoV has been detected and accounted 
for in the error of the estimated Gaussian dispersion, hence al- 
lowing us to properly estimate the uncertainty on the derived 
nuclear NIR luminosity. 



broad distribution, from —0.62 to +0.15 mag dex 1 ) of 12 
nearby bright elliptical galaxies ()Peletier et al.lll990al bh 
ii) -0.35^0 42 m ag dex -1 of a sample of early-type galaxies 
in th e rich cluster AC 118 at z=0.31 ()La Barbera et all 
2002); iii) A(i? - #)/A(logr)=-0.2±0.1 mag dex" 1 for 
the host g alaxies of BL Lac ob jects, thought to be normal 
ellipticals l|Scarpa et alJl200u|) (the H and K host galaxy 
profiles should be similar as dominated by slowly evolv- 
ing red stars). We conclude that based on their colors and 
color gradients the hosts of the studied sample of HEG 
RGs are perfectly consistent with being normal ellipticals. 

3.2.2. IR nuclear measurements accuracy 

Clearly, effects other than nuclear obscuration might cause 
a color change, such as extinction by extended dust struc- 
tures or age/metallicity gradients in the stellar popula- 
tion. As stated above, ~50% of the objects show gradi- 
ents, with only ^20% statistically significant. Such gra- 
dients indeed represent the main limiting factor for the 
accuracy of the color profile method. However, in most of 
the cases showing gradients, the estimates of Lk' based 
on the outer color gradient modeling or assuming a flat 
profile differ by <, 0.2 dex, and never more than a factor 
^2 (see also Appendix 151 for the case of 3CR 403). Only 
3CR 357 and 3CR 379.1 show larger differences, 0.44 and 
0.57 dex, respectively; for the latter one, which has been 
excluded from further analysis, only an upper limit for L^> 
could be estimated. Only 2 (out of 26) objects (3CR 171 2 
and 3CR 284) show complex R — K' profiles. Therefore, 
we can conclude that for most sources the color profile is 
sufficiently well behaved that our rather simple technique 
provides robust estimates of the IR nuclear fluxes, with 
significantly higher accuracy than the surface brightness 
decomposition method. 

Before proceeding any further, it is also crucial to as- 
sess the effect of the uncertainty in the seeing measure- 
ment on the color profile as this is the only free parameter 
of the analysis. In Fig. [5] we compare profiles obtained 
varying the measured FWHM of 3CR 403 by 9% (corre- 
sponding to the variance of the different measurements 
of seeing for this source, slightly larger than that typical 
of the studied sample). The differences between the pro- 
files are only marginal, being less than ± 0.04 in color 
at ff.'5 and completely negligible at 1". As detailed in 
Appendix [BJ this implies an uncertainty on the nuclear 
IR luminosity in the range 0.02-0.2 dex depending on the 
specific source. 



2 The source 3CR 171 is very peculiar, having a very bright 
narrow-line emission r egion co-spatial to the radio emission 
llTadhunter et al1l2000Tl . indicative of jet-cloud interaction; the 
contribution to the HST image fro m line emission is large, if 
not dominant \de Koff et al.ll2000l) . Its observed color profile 
(see Fig. is complex, with a central optical excess likely due 
to the strong line emission contribution (see also Appendix [Q 
for details). We therefore neglect this object in the following 
analysis. 
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Fig. 2. Color profiles of 3CR 403. The filled circles rep- 
resent the profile obtained with the estimated FWHM of 
0'.'89; the large and small empty circles correspond to the 
profiles obtained varying the measured FWHM by 9%, the 
variance of the different measurements of seeing. The dif- 
ference is < 7% in the inner most point, <4% at a radius 
of 0'.'5 and completely negligible at > 1". 
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A crucial test for the robustness of our procedure is 
given by the comparison with the results from direct imag- 
ing for the three sources for which HST-NICMOS obser- 
vations are publicly available, namely 3CR 405, 3CR 456 
and 3CR 460 (fo r deta ils see Appendix 0. For 3CR 405, 
iTadhunter et alj l|l999|) found a nuclear point-like source 
with a luminosity of (4.1±0.2) x 10 41 ergs" 1 at 2.25 /im. 
In the hypothesis of a flat outer R — K 1 profile, we esti- 
mated a nuclear luminosity of (7.3±0.9) x 10 41 ergs" 1 at 4.1. Nuclear properties of narrow and broad lined RGs 



Log L 178 [erg s~ 1 Hz~ 1 ] 

Fig. 3. [OIII] line luminosity, £[oiii]j versus total radio 
luminosity at 178 MHz, Ln%- HEGs are represented by 
filled circles and BLRGs by empty squares. 

by the high spatial resolution optical HST imaging. In 
Table the measured values of the nuclear j/Lk' for the 
HEG sample are listed; the details of the analysis for the 
individual objects can be found in Appendix lO 

4. Results 



2.12 /im, about 0.25 dex larger. However, by modeling 
the small color gradient at radii >1'.'25, we obtained 
i/Lk 1 ^10 41 - 66±0 09 ergs -1 , perfectly consistent with the 
NICMOS result. For 3CR 456 and 3CR 460 NICMOS- 
NIC2 images are available in the public archive at STScI. 
A central point-like source is clearly present in 3CR 456 
with a luminosity of io 43 - 63±009 ergs" 1 at 2.05 /im. After 
modeling the outer color gradient, we estimated vh^ = 
10 43.46±o.o4 ergs -i ; w ithin -2cr from the HST value. In 
the HST image of 3CR 460 only extended emission is visi- 
ble. An upper limit on the IR nuclear luminosity of 10 42 - 82 
ergs" 1 at 2.05 /im has been estimated, consistent with our 
measurement of ^Lk' = 10 42 ' 7 °' 19 ergs" 1 . Considering 
the different methods, aperture and spatial resolution of 
the above observations, the agreement between the HST 
images and the color profile analysis results is quite re- 
markable. Early results from new HST/NICMOS obser- 
vations confirm our measurements for 3CR 79, 3CR 379.1 
and 3CR 403 (M. Chiaberge, priv. comm). 

It thus turns out that the R ~ K' color radial profile 
method returns values of Lk' correct within at most a 
factor of 2, and most likely with an average error of 0.2 dex 
in log, even for very faint nuclei (e.g. 3CR 405). The high 
accuracy basically follows from the information provided 



As pointed out in § 1, in the standard unification picture, 
HEGs are believed to be intrinsically identical to BLRGs, 
and obscuration by a dusty torus is thought to be respon- 
sible for their different nuclear properties. In this scenario 
HEGs should be indistinguishable from BLRGs when com- 
paring isotropic quantities such as the [OIII] line luminos- 
ity, Lroini) ^ n d the total radio luminosity at 178 MHz, 
I/i7g. This indeed appears to be the case: Fig. |3| shows 
Lns versus £rom] for the complete samples of HEGs and 
BLRGs drawn from the 3CR sample within z ^0.3 (from 
CCC02). The Kolmogorov-Smirnofftest and the mean val- 
ues of the two quantities are consistent with the two pop- 
ulations being drawn from the same parent distribution 
(see Tab. |21 for a summary of the statistical comparisons 
between HEGs and BLRGs). 

This unification scenario can now be further tested by 
including the information on the nuclear NIR luminosi- 
ties. In Fig. 21 we have plotted the nuclear optical and 
IR luminosity, Lq and Lk' versus Lns, for both HEGs 
and BLRGs 3 . As already shown by CCC02, BLRGs have 
systematically larger Lq with respect to HEGs: after ex- 



3 We have chosen to use the isotropic Lns instead of L r (as 
done by CCC02) since the latter is likely affected by beaming 
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Table 2. Nuclear luminosities and derived quantities 
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Name 


log(^L K ') a 
ergs -1 




r 


A 

jt-V, torus 




(mag) 


(%) 


(mag) 


(1) 


(2) 


(3) 


(4) 


(5) 



3C 18 


43 66±0 06 


1 87±0 27 


~0 


1.9±0.3 


3C 33 










3C 63 










3C 79 


42.87±0.21 


-0.12±0.87 


5 5+l'5 


60 + ™ 

uu — 29 


3C 98 


41 80+S'nS 


>6 19±0 28 


<0.43 


17 +1 ? 

x 1 - u -1.3 


3C 105 










3C 135 










3C 153 










3C 171 










3C 184.1 


43 16±0 04 


1 67±0 30 


7 + ?2 

' - u —1.7 


19-51?? 


3C 192 










3C 197.1 


42.81±0.07 


0.56±0.38 


~0 


0.6±0.4 


3C 198 


42.57±0.10 


-0.12±0.47 


~0 


-0.1±0.5 


3C 223 


42 43 + °'?5 


>2 34±0 33 


<0.88 


33.2 + 2 ? 


3C 223.1 


43.03±0.04 


>5.08±0.30 


~0 


>5.1±0.3 


3C 234 


44 06±0 04 


1 55±0 30 


6 3+?'° 

U - J -1.6 


21 5 +2 '2 

^ ,i-u — 1.8 


3C 284 


>42.34±0.33 








3C 285 


42.35±0.04 


8.46±0.30 


~f) 


8.5±0.3 


3C 300 


42.54±0.16 


71±0 69 


3 0+?I 


37.7 + 7 2 7 3 


3C 321 


>42.73±0.06 








3C 327 


<41.77 


<2 58±0 29 


<0.22 


44+ S° 


3C 349 


43.02±0.11 


1.15±0.50 


~0 


1.2±0.5 


3C 357 


42.19±0.35 


>2.86±1.41 


<1.0 


29l 10 


3C 379.1 


<42.04 


^1.03±0.28 






3C 381 


42.981™ 


>3.53±0.38 


<1.0 


24.51?;° 


3C 402 


41.59±0.06 


1.71±0.35 






3C 403 


42.27±0.05 


4.17±0.33 


0.88lg;g 


21.9±1.1 


3C 405 


>41.66±0.09 








3C 436 


>42.53±0.36 








3C 452 


42.04l° ;*° 


>2.30±0.47 


<2.5 


24.0l|° 


3C 456 


43.46±0.04 


1.55±0.30 


6.01?;° 


22.01?? 


3C 460 


42.79±0.19 


2.46±0.80 







a log (vL K i): total luminosity at 2.12 /im estimated as in 
6 Nuclear extinction from eq.0 

c Fraction (percentage) of the total nuclear emission scattered toward the line-of-sight (from Fig. |SJ; / ~0 means that the 
emission is consistent with transmitted light both in optical and NIR; 
d Torus extinction (from Fig.|SJ. 



eluding the upper limits (i.e. optically undetected HEGs) 
BLRGs are brighter on average by a factor ^13 (includ- 
ing HEG upper limits increases the difference) . Essentially 
the same result holds for the IR nuclear luminosity (see 
Fig. 01 right panel): the BLRG sub-sample is character- 
ized by larger Lk' 4 with respect to HEGs even though by 
a factor (^ 6) smaller than in the optical band. 

In order to assess whether and to which degree the dif- 
ferent optical/IR luminosities could be ascribed to absorp- 
tion, we estimate the amount of obscuration (Ay) from 
the comparison of the observed ratio ^o/^K'lobs with an 



intrinsic (assumed) value of such ratio, Lo/LK'\ int - More 
specifically: 

A v = —\\og^ -log^- 1, (1) 



where log Lq/Xk/ lint ~-0-42 5 is assum ed equal to that 
typical for RL QSOs l|Elvis et al.lll994|) (and corresponds 
to an optical-to-infrared spectral index of ~0.87). The es- 
timated extinctions are listed in Tabled (column 3): their 
range is < Ay < 4 mag, with only 3 exceptions. 

From these values of Ay an extinction corrected nu- 



(BLRGs show brighter radio cores with respect to HEGs - see clear optical lu minosity Lq. a,, mr r can be inferre d from 
Tab.H. 

4 For BLRGs L K i has been estimated from Lo (from CCC02) 
assuming an optic al-to-infrared spectral index typical of QSOs 
jElvis et alJll994h . 



1 _ — w . ^ w — — ~ — 

A R = 0.751 Ay jCardelli. Clavton fc Mathidll989|) . Even 
after correcting for the estimated nuclear extinction, the 

5 With a typical error of ~0.2 dex jElvis et alJll994h . 
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Table 3. Results of statistical tests: comparison of means 
and standard deviations of different variables for HEGs 
and BLRGs; p is the probability that the two samples are 
drawn from the same parent distribution according to the 
Kolmogorov-Smirnoff test. 



A v [mag] 
4 6 



10 



Sample 


HEG 


BLRG 


P 


log L 178 


33.77±0.48 


33.74±0.50 


0.99 


log £[OIII] 


41.83±0.59 


41.94±0.45 


0.65 


logL r 


30.77±0.75 


31.37±0.42 


0.06 


log Lq 


27.79±0.96* 


28.92±0.57 


6xl0~ 4 * 


log L K ' 


28.57±0.63* 


29.34±0.57 


0.001* 


log I/O, A v corr 


28.38±0.63* 




0.04* 


log I/K', A v corr 


28.80±0.63* 




0.04* 


log vL m 


44.8±0.4 


44.3±0.4 


0.13-0.18 



* only detections for HEGs are considered; in this case, the 
probabilities that the observed log Li7g and log L[oni] distri- 
butions are drawn from the same parent population are p=0.89 
and p=0.57, respectively. 



BLRG optical luminosities systematically exceed those of 
HEGs (see TabEJ). Clearly this still holds for the NIR nu- 
clear luminosity only marginally affected by the relatively 
small Ay (see Tab. |3J). Thus, based on the optical/NIR 
nuclear continuum properties, HEGs and BLRGs do not 
seem to belong to the same parent population even after 
absorption (inferred from the ratio between optical and 
NIR nuclear emission) is taken into account. This seems to 
imply that HEGs could be intrinsically fainter than their 
putative BLRGs counterparts (having similar narrow line 
and extended radio luminosities). 

4.2. Nuclear and emission line properties 

This result appears however to be in contradiction with 
the similarity in the L[om] distributions for HEGs and 
BLRGs, as the same intrinsic nuclear radiation field would 
be expected to maintain similar gas photo-ionization and 
line luminosities (if indeed photoionization by the nuclear 
radiation field is indeed the chief mechanism producing 
L[oui])- 

In such hypothesis, for more quantitative estimates we 
report in Fig. the ratio between vL^' (from Table |2| 
and the ionizing luminosity Li on . L- lon has been estimated 
from the total luminosity in narrow lines Lnlr [inferred 
from £[0111] (Tabled) as Lnt.r = 3(3L[qtt] + 1.5L[q ttt]) 

iRawlings fc SaunderslllQQlj . as 

0.01 6 . 



and Lpm] = 4L[oii]; e.g 



Lion = LnlrC* , adopting a covering factor C 
For BLRGs, log (L ion /pL K ,) ~1.1±0.47. 

When compared with the typical values for BLRGs, 
iion/^Lx' thus represents the 'discrepancy' between the 
nuclear emission as inferred from the lines and the es- 
timated Lk> for HEGs. The representative values for 
BLRGs are indicated in Fig. 5 by the shaded area (ac- 
counting for the dispersion of ~0.47 dex of BLRGs 



/ r 



:?S 7.: : f 



^ 2 



234 m -±'460 
456 W 



O 



3 



4 r 







A 
300 



A sz 381 , 
J367 9 — 

223 



-1 -2 
Log (L /L K .) 



-3 



Note that the precise value of C is ininfluent in the follow- 
ing analysis. 



Fig. 5. Ratio between the ionization luminosity L 1Qn and 
uLk 1 versus the ratio between Lq and Lk'- Symbols as 
in Fig. 0] The arrows represent upper/lower limits. The 
dotted vertic al line indicates t he typical QSO ratio for 
log(Lo/^K') l|Elvis et al.lll994^ . The continuous line rep- 
resents the value Li on /vLK> (estimated for BLRGs) when 
absorption is taken into account [i.e. log {L- mn /vL^i) ~ 
1.1 + 0.048/lv, Ay from eq.d and the shaded area shows 
the location of BLRGs (estimated from the dispersion of 
~0.47 dex around the mean value). The upper axis, la- 
belled Ay, is obtained from the lower one via eq. ^ 



around their mean value shown as continuous line). As 
Lion/vLK' is affected by obscuration [as log (Li on /vLK') ~ 
1.1 + 0.048j4v]j w e plotted this ratio against the observed 
Lo/Lk', which corresponds to Ay (reported on the upper 
axis). Obscuration then 'moves' the position of sources 
in this plane along a line parallel to the continuous one. 
Therefore, if nuclear obscuration is the only cause of the 
differences between HEGs and BLRGs, then HEGs should 
be located within the shaded area, depending on the 
amount of obscuration reflected in their optical-NIR col- 
ors. However, only ~30% of HEGs are consistent with this 
picture. For most sources, instead, the observed L\ m jvLYj 
is much larger than that typical of BLRGs - even when 
extinction is considered - up to factors larger than 20 
(for detections). We conclude that obscuration alone can- 
not account for the differences in the optical/NIR nuclear 
properties of HEGs and BLRGs with respect to the simi- 
larity in their narrow line emission. 

4.3. The role of scattering in HEGs 

An ingredient we have so far neglected is the possible pres- 
ence of a nuclear component due to scattered light. As re- 
ported by CCC02, there are at least two sources (3CR 234 
and 3CR 109) for which spectropolarimetric studies in- 
dicate that their optical nuclear component corresponds 
to a compact scattering region, whose luminosity indeed 
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matches the optical nuclear emission seen in the HST im- 
ages. CCC02 argued that the HST optical cores in most 
HEGs can indeed be produced by scattering. 

The presence of a scattered component radically 
changes the interpretation of our results. Let us examine 
the consequences of this hypothesis in a simple scenario, 
where the observed nuclear luminosity is simply due to di- 
rect light attenuated by a dusty torus characterized by an 
extinction Ak\ torus, and a fraction / of the total nuclear 
light (iK'lint) scattered toward the line-of-sight. E.g. in 
the K' band: 

^| obs =iK'| int 10- a4 ^'.— +/£ K 'L- (2) 
For A r =0.751j4 v and A K ,=0A2A V then: 



L>K' 

and 

Lien 



Lo_ 
L K > 



10 :: 1 ^ ' - + / 



lQ-0.048-Av, 



/ 



obs 



i + / 



;L0-0.048-A v ,toru S _|_ J 



(3) 



(4) 



where log {L /L w ) | int ~ -0.42 and \og(L- ion / vL K ,)\ int ~ 
1.1±0.47 are the typical values for BLRGs. 

The effect of the two added (free) parameters, Ay, torus 
and /, on L-^/iyL^i and Lo/^K' ~~ i-e. the plane defined 
in Fig. - can be graphically seen in Fig. as dotted 
and dot-dashed lines (obtained from eq. Eland E]by fixing 
Av,torus an d letting / run, and vice- versa). As an illustra- 
tive example, let us concentrate on the line representing 
/ = 0.01. For low values of obscuration (A v , torus < 5 
mag) the transmitted light exceeds the scattered fraction 
both in the optical and NIR. A source 'moves' along the 
solid tilted line simply describing the effects of increas- 
ing nuclear obscuration. For larger values of Ay. torus the 
transmitted optical component becomes increasingly neg- 
ligible with respect to the scattered one and Lq is essen- 
tially constant and equal to / Lo|; nt - On the other hand 
Lk' keeps decreasing, i.e. Li ou /vL^i increases, causing the 
turnover toward the bottom left corner of the diagram. 
Only for values of Ay, torus > 40 mag the scattered and 
transmitted if' -band light become comparable: the 'track' 
then halts to the intrinsic Lo /£k' ratio with luminosities 
reduced to a fraction / on the intrinsic ones. 

The properties of all HEGs reported in Fig.Elcan then 
be interpreted with different contributions of transmit- 
ted and scattered light: their location in Fig. translates 
into constraints on / and Ay, torus (listed in Table |21 for 
each source). Only 6 sources, located close to the contin- 
uous line, are consistent with being BLRGs seen through 
a small amount of dust, with Ay. torus < 8.5 mag. The 
majority of HEGs require larger obscuration, namely 15 
< Ay t torus < 50 mag, and their optical cores can only be 
due to scattering. The scattered fraction ranges between 
~ 0.2 and 7%: indeed HEGs with detected optical nucleus 
in the HST images have / > 0.9%, while objects for which 
no optical HST source is seen correspond to upper limits 



in the range 0.2 - 2.5%. Clearly the values of extinction 
derived from the Lq/Lk> ratio (§ 4.1) systematically un- 
derestimate the nuclear absorption. 

Summarizing: within this simple scenario, for the ma- 
jority (~70%) of HEGs both nuclear absorption by a dusty 
torus and scattering of the nuclear light are required to 
unify them with BLRGs. At optical wavelengths, the scat- 
tered component usually dominates (except for the few 
sources which are characterized by a low Ay.torus), while in 
the NIR the nuclear light is generally transmitted through 
the torus, attenuated by a factor typically ~10 but up to 
- 100. 

So far we have neglected the effects of redshift, having 
assumed that all HEGs arc at z=0 in order to avoid intro- 
ducing uncertainties due to the unknown spectral shape. 
In order to quantify these effects on our results, we have 
repeated the above analysis assuming z=0.15, which is the 
median redshift of the BLRG and HEG populations (the 
z distribution of the two population are statistical identi- 
cal, with a Kolmogorov-Smirnov probability p=0.81). At 
z=0.15, the observed optical and NIR luminosities corre- 
spond to rest-frame wavelengths 0.6 and 1.85 /im, respec- 
tively. Consequently, we have re-estimated the rest-frame 
values for Lo/LK>\ int and L lcm /iyLK' | int , and the extinc- 
tion at the appropriate wavelengths. The newly estimated 
values for Ay^orus are smaller by ^25% compared to those 
for z=0 (the rest-frame bands correspond to shorter wave- 
lengths, where the extinction is higher), the fractions / are 
almost unaffected. No trend of ^v, torus with z is present. 
We conclude that our results are not significantly affected 
by the z=0 assumption adopted. 

An important test of our findings is provided by FIR 
information, as the FIR luminosity should be orientation- 
independent. A few HEGs and BLRGs in our sam- 
ple have published FIR (60 um) lu minosities (vL 60 ) 
l|Heckman et al.lll994trHaas et al.ll2004l) . Within the large 
uncertainties and low statistics, the distributions of vLqq 
for HEGs and BLRGs are consistent with being drawn 
from the same population (see Table |3J and the observed 
vLqo/i/Lk' ratio of HEGs is consistent with that of BLRGs 
once the extinction Ay jt0 rus due to the dusty torus is taken 
into account, as shown in Fig. [7| (upper panel). The same 
test ca n be performed using MIR (12 um) luminositie s 
(vL u ) l|Heckman et a.l.l ll994i Isl ^benmorgen e t a,lJ l2004h . 
In Fig. [3 (lower panel) the observed ratio vLyiIvLw is 
plotted instead of vL^jvL^ . Absorption in the MIR has 
been taken into account by using th e extinction la w with 
standard graphite-silicate mixes fe.g. lDraina which 
implies Ayi = 0.036Av. Once the extinction Ay it0 rus is 
taken into account, we find full consistency between HEGs 
and BLRGs (the probability that the vL\i distributions 
of HEGs and BLRGs are drawn from the same parent 
population is p=0A7, when luminosities are corrected for 
extinction, with < vL 12 >=io 44 - 23±0 - 23 for BLRGs and 
< vhxi > = 10 44 - 59±0 - 88 for HEGs). 
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Fig. 7. Upper panel: Ratio between the 60 /jm FIR lu- 
minosity vLqo and vLw (2.12 /mi) versus the torus extinc- 
tion Ay, torus (from Tablc[2J) . Lower panel: Ratio between 
the 12 jim MIR luminosity vL\2 and vL^ versus Ay t torus- 
Symbols as in Fig. [5] The shaded area is the BLRG 
locus: the mean values are vL &0 /vL^ ~0.61±0.74 dex 
and vLyijvIjYi' ~0.67±0.60 dex. The dashed lines in 
both panels represent the BLRG value when absorption 
in the NIR is taken into account [log (vLqq / vLk> ) ~ 
0.61 + 0.048Ay, torus] (upper panel) and when absorp- 
tions both in the NIR and MIR are taken into account 
[log {vLi2/vLw) ~ 0.67+0. 034v4v, torus] (lower panel); the 
two dotted lines delimit the BLRG 'region'. 

4.4. Comparison with previous works 

Three objects (3CR 79, 3CR 223 and 3CR 234) are in com- 
mon with the small sample of detected sources in SWW00. 
The estimated extinctions reported by the above authors 
are consistent with our estimates of Ay, but much smaller 
than t „rn a . lTavlor et al.1 l|l996|) studied a heterogeneous 
sample of RGs and RL QSOs (matched in total radio lu- 
minosity and redshift) using if -band imaging. If the RGs 
and the RL QSOs have the same intrinsic NIR luminos- 
ity (as postulated by the unified scheme), the observed 



difference in their nuclear NIR luminosities could be ac- 
counted for by extinctions Ay ^8-26 mag, with an uncer- 
tainty of ~7 mag, consistent with what found in our work. 
lHaas et alJ l|2004|) have studied the MIR and FIR prop- 
erties of 3CR RGs and quasars, inferring that HEGs are 
consistent with being misaligned BLRGs if their nuclei are 
obscured by extinctions in the MIR o f ~l-2 mag. Using 
the standard extinction law of Drainc (1989), this range 
corresponds to Ay ^27-56, consistent with the range we 
found. 

It is worth also noticing that a search for nuclear 
point-like UV emission in RGs from HST-STIS images 
(NUV-MAMA detector at ~2300 A) has revealed posi- 
tive detection only for one HEG, 3CR 198, implyi ng an 
extinction of Ay=0.08 mag l)Ghiaberge et al.ll2002l) . per- 
fectly consistent with the value estimated in this work 
(-0.12±0.47 mag). None of the other HEGs targeted 
(3CR285 3CR 321 and 3CR 452) show point-like UV 
emission ((Ghiaberge et al.(l2002l) . For the only object with 
an optical detection (3CR 285), as suming ^4tjv=2.26 Ay 
l(Cardelli. Glavton fc Mat hislll989 |) 1 an o ptical-UV spec- 
tral index 1.0±0.3 llChiaberge et alJl2002l consistent with 
that of RL QSOs. lElvis et all ll994jT and a UV flux 
limit o f ~1.9xl0~ 19 ergs~ 1 cm~ 2 A (the lowest detec- 
tion in IChiaberge et al.ll2002j) . we estimate an extinction 
Ay ^3.5 mag, consistent with our previous finding. 

4.5. Obscuration vs orientation ? 

A key ingredient that should be added into this analysis is 
the source orientation. In the simplest unified scheme, the 
system is supposed to be axisymmetric, with the jet direc- 
tion corresponding to the symmetry axis. The dusty torus, 
perpendicular to it, would hide the central quasar for an- 
gles with the linc-of-sight larger than the torus opening 
angle. 

In principle an estimate of the inclination could be 
derived from the observed correlation between the radio 
core luminosity a t 5 GHz and the radio total luminos- 
ity at 408 MHz l|Giovannini et all l200lh . e.g. by com- 
puting the normalized core power Pcn, defined as the 
source core power normalized to that corresponding to 
an average observing a ngle of 60° for a given total power 
l|de Ruiter et alJll990l) . Within the simplest scheme, one 
would then expect an anti-correlation between Pcn and 
Ay, torus- Although there seems to be a broad trend for 
HEGs with lower Pcn to have higher extinction (the trend 
is stronger if BLRGs are also considered), the large un- 
certainties and low statistics do not allow to draw any 
significant conclusion. 

5. Discussion 

From the above results, the following picture emerges for 
HEGs. On one hand, based on their extended and isotropic 
properties, HEGs are consistent with being drawn from 
the same parent population of BLRGs. In particular, 
HEGs and BLRGs have the same distributions of Ln 8 
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and L[oiii]- If the narrow-line emission arises from photo- 
ionization by the nuclear radiation field, the similarity in 
£[om] implies that the intrinsic nuclear luminosities of 
these populations should also be similar. On the other 
hand, based on their nuclear optical and NIR luminosi- 
ties, HEGs and BLRGs do not appear to be taken from 
the same parent population, with HEGs being signifi- 
cantly underluminous. Except for a handful of sources 
(3CR 18, 3CR 197.1, 3CR 198, 3CR 223.1, 3CR 285 and 
3CR 349), obscuration alone cannot account for such dif- 
ferences. Moreover, for the great majority of HEGs the 
observed nuclear NIR luminosity is insufficient to account 
for the narrow-line emission. 

However, these findings can be reconciled with the uni- 
fication scenario for HEGs and BLRGs if a scattered com- 
ponent is present. Indeed CCC02 have suggested that the 
optical nuclear emission in HEGs can actually be identi- 
fied with a compact scattering region. It is thus possible 
that for the majority of HEGs the observed nuclear optical 
emission is dominated by the fraction of the nuclear light 
scattered toward the line-of-sight. At NIR wavelengths, in- 
stead, both a transmitted (direct) and the scattered com- 
ponent contribute to the observed luminosity. In this (sim- 
ple) interpretation the fraction of scattered nuclear light 
ranges from a few tenths of to a few percent, while the 
torus extinction Ay, torus spans ~15 to 50 mag, with only 
a few exceptions. 

This scenario is consistent with independent findings 
on the best studied object of our sample: 3CR 234. 
3CR 234 is characterized by a very high degree of polar- 
ization (~17% after starlight subtractio n) of both the op - 
tical continuum and the broad Ha line l(Antonuctilll984h . 
which can be explained only with polarization induced 
by scattering; also the -fT-band c ontinuum is polari zed, 
at t he lower level of ~ 4.6±0.4% l|Sitko fc ZhiJll99li see 
also lYoung et al.lll998h . The flux and polarization prop- 
ert ies at optical and N IR wavelengths have been modeled 
bv lYoung et al.l £998), who concluded that the observed 
light comprises a scattered component and a dichroic view 
through the torus which becomes important only at the 
longer wavelengths. This is qualitatively consistent with 
our results, although the nuclear extinction they estimated 
(~60 mag) is larger than ours (22±2 mag). 

Among the sample of HEGs we found 6 sources requir- 
ing low obscuration. Four of them (3CR 18, 3CR 197.1, 
3CR 198, 3CR 349) are consistent with Ay, torus <2 mag, 
i.e. essentially no obscuration, apparently in contrast 
with the lack of broad lines in their optical spectra. 
Interestingly, results obtained during the preparation of 
this work indicate that 3CR 18 and 3CR 197.1 are actu- 
ally BLRGs, misclassified as HEGs. A high S/N VLT spec- 
trum revealed a prominent broad (FWHM^4500 km s _1 ) 
Ha line in 3CR 18 and a broad Ha line is clearly visible 
in the SDSS spectrum of 3CR 197.1. The hypothesis that 
misclassification affects also 3CR 198 and 3CR 349 ap- 
pears more likely than the alternative explanation that 
they actually lack a BLR. Conversely, 3CR 223.1 and 
3CR 285 show sufficient obscuration to prevent the de- 



tection of broad lines in the optical (e.g. the Ha would 
be reduced by a factor > 40), but should be easily de- 
tected in the NIR. The presence of these low-extinction 
objects among the sample is particularly interesting with 
respect to whether there is continuity between HEGs and 
BLRGs - i.e. if the torus optical depth varies smoothly or 
if a minimum threshold exists. Higher quality optical and 
NIR spectra are clearly needed to settle this issue. 

The population of HEGs can thus be characterized by 
a distribution in the amount of absorbing dusty material, 
although it is not possible from the current data to deter- 
mine whether this is closely associated to the orientation 
(i.e. the vertical structure of the torus) or simply to differ- 
ent sources having different amount of dust. In any case, 
for a handful of HEGs the direct view of the central engine 
is possible both at optical and NIR wavelengths (might be 
because the line-of-sight intercepts only the edges of the 
torus). At higher levels of extinction the nuclear optical 
light from the central engine is absorbed. The observed 
point-like optical emission is due to a compact scattering 
region located 'outside' the obscuring torus. At NIR wave- 
lengths, instead, the central nucleus is still directly visible, 
and both the direct /transmitted light passing through the 
torus and the scattered components contribute to the ob- 
served NIR luminosity. This is the situation for the major- 
ity of HEGs. The possibility that even higher extinction 
is present is some sources remains open (e.g. 3C 379.1 ?), 
although their number is extremely limited according to 
our statistics (note that for very high extinction also ra- 
diation re-processed by the dust could contribute to the 
observed emission). 

Let us finally compare the ranges of values of the 
scattering fraction and obscuration we inferred with 
respect to those estimated in the literature for ob- 
scured AGN (also non RL, namely Seyfert 2 galax- 
ies). The scattered fraction in the latter ones spans 
from less than ~0 1% (e.g. NGC 1068 iMiller fc GoodrichI 
Il99lt 3CR 321. lYoimg et al.l Il996al 1996b) un to a 
few % (e.g. Idi Se rego Alighieri. Cimatti fc Fosbur\1ll994l: 
iKishimoto et al.l2002aHMaiolino. Marconi fc Olival2000|V . 
A rather conservative upper limit of / < 10% has been 
estima ted for the Seyfert 2 Mrk 477 l|Kishimoto et alJ 
2002b). Our values of / are thus broadly consistent with 
those reported for Seyfert 2, although only for a very 
few sources this could be truly determined. Also the in- 
ferred values for the torus extinction appear to be con- 
sistent with those observationally and theoretically es- 
tima ted for Seyfert 2. For example, some t orus models 
(e.g. iGranato. Danese fc Franceschinil Il997l iFadda et alJ 
1998) require - in order to reproduce their IR spectra 
- moderate-to- high visual extinctions (Ay = 10-80 mag). 
The d etection of NIR broad emission li nes in some Seyfert 
2 (e.g. lVeilleux. Goodrich fc Hm1ll997h also indicates that 
at these wavelengths the nuclei can be directly seen, with 
estimated Ay from ~8 mag up to more than 68 mag. 



12 



Marchesini et al.: Probing obscuration in radio-galaxies 



6. Conclusions 

We have presented i-T'-band observations of all 3CR HEGs 
(FR II) at z <0.3, providing for the first time homoge- 
neous NIR imaging of a complete sub-sample of RGs. After 
showing that the technique based on the surface bright- 
ness decomposition to measure central point-like sources is 
affected by large uncertainties for the objects in the stud- 
ied sample, we have introduced a more accurate method 
exploiting the R — K' color profiles to measure the NIR 
luminosity of a central point-like source. In all but two 
HEGs, a substantial if'-band excess is found, most likely 
associated to the nuclear emission. 

The NIR (2.12 /im) nuclear luminosities have been 
measured for the sources, and used to test the unifica- 
tion scheme for RL AGNs. The picture emerging from 
this work is that for the majority of the studied HEGs, 
obscuration alone cannot account for the different nuclear 
properties of HEGs and BLRGs, and also scattering of the 
(optically) hidden nuclear light toward the line-of-sight 
should play a role. For ~70% of the HEGs, the observed 
point-like optical emission is dominated by such scattered 
component (the direct one being totally absorbed), while 
in the -fT'-band both transmitted light and the scattered 
one contribute to the observed nuclear luminosity. The es- 
timated fraction of scattered light ranges from a few tenths 
of to a few percent of the hidden AGN light, while the 
torus extinction is in the range 15< Ay t torus <50 mag, 
with only a few cases with lower obscuration. Interestingly, 
two of the four HEGs consistent with low levels of obscu- 
ration (Ay, torus < 2 mag) turned out to be mis-classified 
objects in which a prominent BLR emerges in high quality 
spectra. This hints to the possibility that the remaining 
two low-obscuration objects are also mis-classified BLRGs. 

To further test this scenario, deep NIR spectropo- 
larimetric observations are fundamental for two reasons. 
Firstly, they should allow to detect hidden broad lines 
even when affected by a certain amount of extinction. 
Secondly, a measurement of their polarization can give 
an estimate of the scattering fraction. Another important 
piece of information could be provided by polarimetric, 
multiwavelength observations, as they would allow to dis- 
criminate between electron and dichroic scattering. HST- 
NICMOS observations will be clearly of key importance 
to determine the presence and measure the luminosity 
of nuclear point-like sources. Finally, X-ray spectroscopy 
can provide estimates of the nuclear extinction: assum- 
ing a standard gas-to-dust r atio (Ay = 5 x lO -22 ^, 
iBohlin. Savage fe Drakell978|) . we expect the HEGs in our 
sample to be Compton-thin, with column densities in the 
range A H =0.3-1 x 10 23 cm" 2 . 
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Appendix A: Surface brightness decomposition 
(GALFIT) 

We have applied the 'standard' surface brightness decom- 
position technique to detect and measure point sources 
at the centre of galaxies b y using the two -dimensional 
fitting algorithm GALFIT l)Peng et al.ll2002(l . Th e host 
galaxy surface brightness profile was modeled with a lSersid 
( 1968) profile; the nuclear point-like component was mod- 
eled with a Moffat function, the parameters of which have 
been constrained using the stars present in the FoV. 

In Fig. IA.1I the results of the GALFIT analysis are 
shown for 3CR 456, 3CR 234 and 3CR 18. These three 
objects have the largest HST optical nuclear luminosi- 
ties, vLq >10 43 ergs -1 . Moreover, 3CR 456 (as detailed 
in Appendix has been observed with HST-NICMOS, 
which clearly reveals the presence of a central point-like 
source with luminosity ^io 43 - 63±0 - 09 ergs -1 at 2.05 [m\ 
(F205W filter). The GALFIT best fit solution for 3CR 456 
returns a magnitude for the central point-like source 
mK'=14.90±0.07, which translates into a nuclear lumi- 
nosity of io 44 06±0 03 ergs -1 at 2.12 fim. This is -0.6 dex 
larger than the value from the HST image. For 3CR 234, 
the GALFIT analysis gives mK'=13.86±0.07, implying a 
luminosity of io 44 - 25±0 03 ergs -1 . This value is consistent 
wit h the one measured by SWW00 of — 10 44 32 ergs -1 and 
bv iTavlor et alJ (|l996|) of — 10 44 2 ergs -1 from ground- 
based if -band observations. As detailed in Appendix El 
the value measured with GALFIT is also consistent with 
that obtained with the color profile technique. Finally for 
3CR 18, GALFIT detects a central point-like source with 
magnitude mK' = 14.97±0.07, i.e. a nuclear luminosity of 
_ 10 43.83±o.O3 ergs -i. As for 3CR 234, this value is con- 
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Fig. A.l. GALFIT imaging decomposition for 3CR 456, 
3CR 234 and 3CR 18. The left-hand panel is the observed 
field; the central one is the model fitted by GALFIT; the 
right-hand panel is the residual image. The circles indicate 
the fitted sources. 



sistent within a factor ~1.5 with the nuclear luminosity 
measured via the color profile. 

In Fig. [01 the results from the GALFIT analy- 
sis for 3CR 402 and 3CR 405 are shown. These ob- 
jects have been chosen because the former is the near- 
est source (z= 0.025), while 3CR 405 has been observed 
with NICMOS l|Tadhunter et al.lll999|) . For both sources, 
we performed the GALFIT analysis without and with 
a central point-like source (case A and case B, respec- 
tively). For 3CR 405, the NICMOS image (F222M filter) 
clearly shows a point-like source with a nuclear luminos- 
ity of (4.1±0.2)xl0 41 ergs" 1 at 2.25 urn l|Tadhunter et alJ 
1999). The GALFIT analysis performed without adding 
a point-like source is identical to that with a central 
source with nuclear luminosity a factor ^2.5 larger than 
the NICMOS one: based on the goodness of the fit, it is 
impossible to detect the presence of a nuclear source in 
the ground-based image. Moreover, the best-fit brightness 
of the point-like source measured with GALFIT depends 
on the initial guessed value, with differences as large as 
3.2 mag between best-fit models with the same goodness. 
A similar result is obtained for 3CR 402, for which the 
quality of the best-fit models with and without a central 
point-like source is the same. 

We conclude that the standard technique to detect and 
measure point-like sources at the centre of galaxies (per- 
formed with GALFIT) returns stable and consistent re- 
sults only when the point-like source dominates over the 
host galaxy (as for 3CR 456 and 3CR 234). This is not the 
case for the great majority of the objects in the studied 
sample. As shown for 3CR 405 and 3CR 402, the standard 
technique cannot be used for faint point-like sources, even 
for the nearest objects, as the GALFIT analysis becomes 
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Fig. A. 2. GALFIT imaging decomposition for 3CR 405 
and 3CR 402. Panels as in Fig. IA.ll Case A does not in- 
clude in the modeling analysis a central point-like source, 
while case B does. 

unable to detect such component and the resulting lumi- 
nosities become strongly dependent on the initial values. 

Appendix B: The R — K' color profile technique 

The R — K' color profile technique is here described in de- 
tail. We consider two representative sources. For 3CR 403, 
a fair example of the 'good quality' sample, we illustrate 
the whole technique, i.e.: estimate the nuclear IR lumi- 
nosity in the hypothesis of a 'well-behaved' color profile, 
namely a constant outer (galaxy) R — K' profile; show how 
the uncertainty on the seeing measurement affects the fi- 
nal result; repeat the analysis in the more realistic case in 
which the host galaxy shows color gradients. We then con- 
sider 3CR 300, representative of the 'poor quality' sample, 
having a rather noisy color profile. 

B.l. 3CR 403 

As explained in § [3J the basic idea is that the presence 
of an obscured nucleus will reveal itself with an increase 
of the IR nuclear flux with respect to what is seen in the 
optical images. In order to compare the NIR and optical 
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Fig. B.l. R — K' color radial profile (left-hand panel) and 
aperture flux radial profile (right-hand panel) for 3CR 403. 
The dashed vertical line represents the dispersion (FWHM 
divided by 2.354) of the point-spread function of the im- 
age. The error bars indicate la errors. 

images, it is necessary to produce a synthetic HST image 
that matches both the seeing conditions and the pixel size 
of each IR image. We fitted a Gaussian to all objects in 
the HST FoV which are clearly stellar and measured their 
FWHM. The adopted FWHM is their median value and 
its uncertainty was estimated from their dispersion. For 
3CR 403, FWHM=0'.'89±0.08, and therefore the error on 
the seeing is ~9% - slightly larger than that typical for 
the studied sample. Each HST image has been then con- 
volved with a Gaussian with the appropriate FWHM and 
interpolated to match the pixel size of the TNG image. On 
the 'matched' HST image and on the reduced TNG one 
we performed photometry on annuli of increasing radius 
out to 3"and an average color (and its associated error) 
was derived for each annulus, providing us with the radial 
color dependence of the source. 

The R — K' color radial profile, shown in Fig lB.ll (left- 
hand panel), is obtained from: 

R- K' = 2.5 log 1^^-) +4.025, (B.l) 

y FR jan j 

where Fk'.an and -pR, an are the annulus fluxes in the TNG 
and HST image, respectively; the corresponding error is 
estimated with the standard error propagation. 

To estimate the nuclear monochromatic FT'-band lumi- 
nosity Lk' , the FT'-band aperture flux ) radial profile 
is also needed and estimated from the counts within the 
circular aperture (see Fig. IB. II right-hand panel). 

For a 'well-behaved' R — K' profile, Lk' can be then 
estimated with the following steps: 

1. the reference host galaxy 'color', Fk> / Fri| gal , is mea- 
sured in the outer part of the profile (at r=2'.'7b for 
3CR 403); 

2. the aperture-corrected nuclear IR excess luminosity 
Lk',xs is obtained by multiplying the measured nuclear 
IR excess [defined as I - [(F K ,/F R \ g J/(F K ,/F R )}, 
where Fk' / Fr is the nuclear color] with the corre- 
sponding Fk' • 

3. the HST optical nuclear component Lq is properly 
scaled with the host galaxy color FK'/FR.| ga i- 

4. Lk> is finally estimated as Fk' = ^K'. scaled + Fk',xs- 
The error on Lk' has been inferred from the standard 



Table B.l. IR luminosities for 3CR 403: fixed galaxy 
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modeled galax 


color 








FWHM 


log-t/ K ',xs 


log -^K', scaled 


log vL w 


(") 


(ergs" 1 Hz- 1 ) 


(ergs" 1 Hz" 1 ) 


(ergs -1 ) 


0'.'89 


28.05±0.05 


27.26±0.07 


42.27±0.05 


O'.'89+la- 


28.19±0.04 


27.25±0.07 


42.39±0.04 


0'.' 89- Ict 


27.84±0.09 


27.27±0.07 


42.09±0.08 



error propagation, taking into account also the error 
on L (-10-20%). 

The above procedure has been repeated twice to take 
into account the uncertainty on the seeing determination: 
the HST image has been convolved with a Gaussian with 
a FWHM increased and decreased by its la uncertainty 
and the whole analysis repeated. In Table lB~T1 the derived 
luminosities are listed. The difference in the final Lk' for 
the different seeing is within the la error and of order 
-0.07 dex. 

As shown in Fig. IB.2| f upper panel), there is a diffuse 
IR excess at radii ^l'/25 due to a small color gradient, ne- 
glected so far. Such gradient has been fitted with a linear 
model R — ^'| mode i = A + Br at radii <T'.'25, and extrap- 
olated at smaller radii [the profile and the fitted model are 
shown in Fig. IB.2| f middle panel)]. This best-fit model is 
then used as the host galaxy reference Fk' /FR,| gal at each 
radius (the associated error has been taken as the scat- 
ter around the best-fit) to derive the effective IR excess 
(Fig. IB. 21 lower panel) . In Table IB. 21 the estimated lumi- 
nosities are listed, also showing the effects of the seeing 
uncertainty, which is now of order ~0.15 dex. 

The color gradients indeed represent the main limiting 
factor for the accuracy of this method. However, in most 
of the cases showing gradients, the estimates of Fk' based 
on a flat profile or a model of the gradient differ by 
0.2 dex, and never more than a factor —2 (only 3CR 357 
and 3CR 379.1 show larger differences, 0.44 and 0.57 dex 
respectively). 

B.2. 3CR 300 

In Fig. IB. 31 (left-hand panel), the radial color profile for 
3CR 300 has been plotted (showing also the effects of 
the seeing measurement uncertainty, —2% for this source). 
The profile is quite noisy, and while it is consistent with 
being constant, a color gradient seems present. For this 
reason, we modeled the color profile as in § IB. II (the color 
profile and the fitted model arc shown in Fig. IB. 41 left- 
hand panel). 



16 



Marchesini et al.: Probing obscuration in radio-galaxies 



-0 
-0 

-0 

uT -0 

-0 

-0 
-0 



: : i 
O 



1 .0 

0.8 
0.6 
0.4 
0.2 
0.0 


20 

25 

30 

35 

40 

45 
50 


1 .0 
0.8 
0.6 
0.4 
0.2 
0.0 














- 3.2 


< 


-16.6 
-16. a 




\ 

„ I 


- 3.0 


E 


-17.0 




I 

I 

I 


k 










-2.8 J, 


cp 

QJ 


-17.2 
-17.4 




I 

I 


- 2.6 


Lc" 

!: 


-17.6 







0.0 C.5 1.0 1.5 2.0 2.5 3.0 
R [arcsec] 



0.0 0.5 1.0 1.5 2.0 2.5 3.( 

R foresee] 



JJJLU-J--H4 



0.5 1.0 1.5 2.0 2.5 3.0 
R [arcsecl 



Fig. B.3. Color profile (left-hand panel) and aperture flux 
profile (right-hand panel) for 3CR 300. Symbols as in 
Figs. H and EU 
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Fig. B.4. Color profile with over plotted the fitted model 
(left-hand panel) and the derived IR excess (right-hand 
panel) for 3CR 300. The best-fit line in the left-hand panel 
is R— iT| modcl =(3.11±0.06)-(0.20±0.10)r, and the sco- 
ter of the points around the best-fit line is 0.02 mag. 







Table B.3. IR luminosities for 3CR 300: 


modeled galaxy 






color 










FWHM 

(") 


log L K > , xs 
(ergs" 1 Hz" 1 ) 


log Lk', scaled 

(ergs" 1 Rz- L ) 


log vL K > 
(ergs -1 ) 






0'.'68 


27.151^ 


28.36±0.07 


42.54±0.16 






(y.'68+lcr 


2726 +CK78 


28.37±0.07 


42.55±0.17 






O'.'68-lcr 


27 13 +085 


28.35±0.07 


42.53±0.16 


i 


i 


Appendix 


C: Individual galaxies 
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Fig. B.2. Upper panel: the radial profile of the IR ex- 
cess for 3CR 403 in the hypothesis of a constant host 
galaxy color; at radii larger than l'/25 a small diffuse IR 
excess indicates that the host color profile is not perfectly 
flat. Middle panel: the observed color radial profile with 
over plotted the best-fit model to the color gradient at 
radii larger than l'/25. Lower panel: as the upper panel 
after modeling the small (but significant) color gradient: 
only the nuclear IR excess is present (while the diffuse 
one is no more present). The vertical dashed line is as in 
Fig. IB. II the dot-dashed line represents the zero IR excess 
line. Error bars are la errors. 



In Table IB73I the estimated luminosities for 3CR 300 
are reported. After modeling the color gradient, no signif- 
icant IR excess is found (see Fig. IB .41 right-hand panel) 
and therefore L K > — ^K', scaled- 



Here we report the analysis of the R — K' radial pro- 
file and the results used in the estimate of the nuclear 
IR luminosity for the individual sources. In Tab. IC.ll we 
listed the seeing, the IR excess, the reference host galaxy 
color R — K' | gal , the parameters of the linear model fit- 
ted to the outer color profile to account for possible gra- 
dients, and short comments. We also comment on indi- 
vidual sources when needed and compare our results with 
the previous ones - based on the standard technique - by 
SWWOO and lTavlor et~all Jl99fih . 

3CR 079: SWWOO measured vL K = 10 43 - 66 ergs" 1 , a fac- 
tor of ~4-6 larger than our measur ement, and a facto r 
3-5 fainter than the one reported in iTavlor et alJ^199fi|> . 
The discrepancy between these measurements is not eas- 
ily understandable. While the observed flat R — K 1 profile 
excludes a possible underestimate of Lk', the standard 
technique requires a more uncertain modeling of the host 
galaxy, possibly causing to underestimate the host contri- 
bution. 

3CR 098: SWWOO measured uL K <10 4128 ergs" 1 , while 
we obtain a luminosity a factor of ~3 larger. 
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Table C.l. Details of the color profile technique for the individual sources 
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Source 


Seeing 


IRE 


gal 


A 


B 


Comments 


Name 


(") 


(10 28 erg/s/Hz) 


(mag) 


(mag) 


imag_\ 
\ " 1 




3C 


18 


1.43±0.20 


11.3±2.9 


3.51±0.19 


3.19±1.86 


0.22±1.12 




3C 


33 


0.87±0.11 












3C 


79 


1.68±0.03 


2.40±2.26 


3.09±0.02 


2.89±0.29 


0.17±0.16 


if=14.79±0.08 a ; if=14.79±0.06 (SWW00) 


3C 


98 


1.71±0.06 


0.43±0.03 


2.93±0.01 


2.89±0.02 


0.03±0.01 


if=14.71±0.08; if =12.60±0.03 (SWW00) 


3C 


105 


1.73±0.08 












3C 


153 


2.23±0.07 












3C 


171 


1.02±0.04 


-9.39±2.83 


2.30±0.14 






if=15.69±0.12; if =15.71±0.04 (SWW00) 


3C 


184.1 


1.26±0.06 


3.73±0.31 


3.16±0.05 


3.23±0.16 


-0.05±0.09 


well-behaved color profile 


3C 


197.1 


0.87±0.07 


0.49±0.29 


2.87±0.02 


2.87±0.06 


0.01±0.05 


well-behaved color profile 


3C 


198 


0.97±0.03 


-0.14±0.40 


2.62±0.08 






small nuclear IR deficit 


3C 


223 


1.36±0.06 


1.35±0.21 


2.86±0.03 


2.72±0.17 


0.13±0.11 


if=14.54±0.08; if =14.58±0.05 (SWW00) 


3C 


223.1 


0.73±0.11 


7.00±0.71 


2.87±0.03 






well-behaved color profile 


3C 


234 


0.90±0.05 


56.1±5.4 


2.31±0.13 






if=13.48±0.08; if =13.49±0.05 (SWW00) 


3C 


284 


0.63±0.07 


1.55±1.19 


2.99±0.06 






complex profile with a small bump 


3C 


285 


0.89±0.04 


1.55±0.13 




3.13±0.02 


-0.12±0.01 


largest outer color gradient 


3C 


300 


0.68±0.02 


0.14±0.27 


2.79±0.24 


3.11±0.06 


-0.20±0.10 


noisy color profile 


3C 


321 


0.92±0.13 


3.68±0.17 


2.60±0.06 


2.63±0.17 


-0.01±0.08 


well-behaved color profile 


3C 


327 


0.61±0.06 


<0.28 


3.06±0.04 






flat color profile at r ^1" 


3C 


349 


0.81±0.03 


0.55±1.49 


3.15±0.10 


3.38±0.05 


-0.22±0.07 


noisy color profile 


3C 


357 


0.71±0.10 


0.76±0.84 


2.97±0.12 


3.16±0.03 


-0.09±0.03 


noisy color profile 


3C 


379.1 


0.67±0.03 


<J2.16±6.57 


3.22±0.21 


3.36±0.09 


-0.12±0.10 


very noisy color profile 


3C 


381 


0.80±0.11 


5.40±1.03 


3.12±0.15 


3.10±0.40 


0.01±0.27 


well-behaved color profile 


3C 


402 


0.77±0.08 


0.15±0.03 


2.82±0.02 






well-behaved color profile 


3C 


403 


0.89±0.08 


1.13±0.14 


3.01±0.05 


3.15±0.04 


-0.05±0.02 


well-behaved color profile 


3C 


405 


0.95±0.04 


0.32±0.07 


3.29±0.03 


3.46±0.02 


-0.06±0.01 


well-behaved color profile 


3C 


436 


0.72±0.20 


2.37±1.94 


3.33±0.23 


3.61±0.10 


-0.20±0.09 


noisy profile 


3C 


452 


0.96±0.04 


0.46±0.17 


3.04±0.04 


3.14±0.07 


-0.05±0.04 


well-behaved color profile 


3C 


456 


1.00±0.07 


10.8±1.0 


2.92±0.07 


2.64±0.41 


0.28±0.36 


if=14.64±0.06; if =14.63±0.10 (dV98) 


3C 


460 


0.88±0.10 


^2.20±6.49 


3.31±0.17 


3.46±0.31 


-0.04±0.24 


if=14.69±0.08; if =14.71±0.10 (dV98) 


IRE 


is the nuclear aperture-corrected infrared (2.12 /im) 


excess Lk' 




is the host galaxy color; A and B are the 



parameters of the linear model R — K'\ mode] = A + Br fitted to the observed outer color profiles to account for possible color 
gradients [when not specified, the outer profile is consistent with being flat (3CR 198, 3CR 223.1, 3CR 234, 3CR 402) or 
complex (3CR 171, 3CR 284)]; a if is the if —band phot ometry within a 3"-aperture for the objects from SWW00, and within 
a /.^-aperture for the sources from Ide Vries et all <|l998h (dV98). 



3CR 171: A very significant gradient is found, from 
r=l'.'75 to the innermost point corresponding to a nuclear 
IR deficit of (9.39±2.83)xl0 28 ergs" 1 Hz" 1 at r=0'/5. 
This might be due to an optical excess caused by the 
strong line emission presen t in the inner A" (Ha flux of 
~2.12xl0" 14 ergs -1 cm" 2 : lTadhunter et al.ll2000^ . 

3CR 223: SWW00 measured vL K ~ io 43 - 50 ergs" 1 , a 
factor of larger than ours. 

3CR 234: Our 3" aperture photometry is about 0.4 mag 
larger than the value reported in SWW00. However, 
using their best fit parameters for the host+nucleus 
decomposition, we estimate a 3" aperture photometry 
of K =13.49±0.05, perfectly consistent with our value. 
SWW00 measured i/L K =10 44 ' 32 er gs" 1 , a facto r ~1.8 
larger than our measurement, while rTavlor et aTl l)l996|) 
measured ^Lk=10 44 ' 2 ergs" 1 . These two values are con- 
sistent within the errors, and this is not surprising as this 
source is dominated by the nuclear emission, and are also 
consistent with our value. 



3CR 327: The R — K' profile is characterized by a gradi- 
ent at radii >1" and a constant R — K' color toward the 
center (R — K '=3.06±0.04 applying a weighted mean of 
the innermost 8 points). The inner plateau implies a la 
upper limit of the nuclear IRE of 2.77xl0 27 ergs" 1 Hz" 1 . 
3CR 349: The color profile is quite noisy, and the sky 
dominates the profile at radii >l'/75. Assuming a flat outer 
(r ^0"75) color profile, we estimated vL K , =io 43 - 09 ±°- 10 
ergs" 1 . Although very noisy, the color profile shows a gra- 
dient in the outer part; we therefore studied the extreme 
case in which no IRE is present by modeling the color 
profile. The resulting vL^j = is io 43 02±011 ergs -1 , con- 
sistent with the previous estimate. 

3CR 357: Assuming a flat outer color profile, we esti- 
mated vLy^' = io 4263±0 - 21 ergs" 1 ; modeling the outer 
(r >0"4) galaxy profile, we obtained vL K , = io 4219±0 - 34 
ergs -1 , 0.44 dex smaller than the former estimate. 

This object and 3CR 379.1 are the worst ones, with vL^* 
badly constrained (but still within an uncertainty of a fac- 
tor ~3). In the HST image, extended dust lanes are clearly 
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visible within the galaxy, down to the very inner region. 
In order to test the sensitivity of our analysis to extended 
dust lanes, we modeled the galaxy in the HST image with 
an elliptical model (masking out all the visible dust lanes). 
Then we repeated our analysis using the 'dust-free' mod- 
eled galaxy profile: the two R — K' profiles are very simi- 
lar (differences well within the error bars). Therefore, we 
conclude that our analysis is not significantly sensitive to 
extended dust lanes. 

3CR 405: The observed R — K' profile is very well- 
behaved, being flat at radii >2" with variation <0.01 
mag. The nuclear IRE measured is (5.13±0.64) x 10 27 
ergs -1 Hz -1 . 3CR 405 (a.k.a. Cv gnus A) has been pre- 
viously observed with NICMOS feadhunter et al.lll999h 
(F222 M filter centered at 2.25 urn). iTadhunter et al.1 
(1999) report the presence of a nuclear point source (un- 
resolved, with FWHM <0'/21) which dominates the emis- 
sion at 2.25 ^m with a flux F 2 .25 Atro =(4.9±1.0) x 10 -28 
erg s -1 cm -2 Hz -1 . This translates into a luminosity of 
(4.1±0.2)xl0 41 ergs -1 at 2.25 /im, about a factor 
~1.8 smaller than the value reported in this work of 
(7.3±0.9) x 10 41 ergs" 1 at 2.12 /im within an aperture of 
~0'/75. Our larger value is due to the fact that we did ne- 
glect the color gradient at r iH'.'25: if this gradient is fitted 
with a linear model we obtain i/Lk' = (4.57±0.94) x 10 41 , 
perfectly consistent with the NICMOS value. 
3CR 456: The F205W filter image from the STScI public 
archive clearly shows the presence of a nuclear point-like 
source with i/Lj,=10 43 - 63±0 09 ergs -1 at 2.05 /im (the er- 
ror, about 20%, is dominated by the uncertainty on the 
subtracted background) consistent at the 1.8a level with 
our measurement. 

3CR 460: 3CR 460 has been observed with NICMOS 
(F205W filter). From the NIC2 image in the public archive 
the emission is extended, and there is no clear sign of a nu- 
clear point-like source. We set as an upper limit the light 
excess of the central circular aperture of radius 0'.'17, ob- 
taining vL v <io 42 - 82±0 04 ergs -1 at 2.05 /im, consistent 
with our measurement of io 42 - 79±019 ergs -1 at 2.12 /xm. 



Appendix D: Optical and NIR images 
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Fig. 1. Color radial profiles of the 26 HEG RGs [R - K' = 2.5 log (F K > /F R ) + 4.025]. The error bars represent the la 
error on the derived color, while the dashed vertical line represents the dispersion (FWHM divided by 2.354) of the 
PSF of the image. 
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31 32 33 34 35 36 32 33 34 35 36 

Log L 178 [erg s~ 1 Hz~ 1 ] Log L 178 [erg s _1 Hz~ 1 ] 

Fig. 4. Left panel: Optical core luminosity Lq vs Lns- Right panel: NIR nuclear luminosity Lk' vs Lyj%. BLRGs 
are represented by light filled diamonds; the different symbols to represent HEGs correspond to the different quality 
of the R — K' color profile: filled circles for good quality, triangles for worst quality and empty circles for intermediate 
one. Typical error bars for the three quality levels are also shown; arrows are upper limits and the short vertical lines 
associated to the few objects in the right-hand panel indicate the value of Lk' if their Lo where equal to zero (since 
the latter are only upper limits). Squares represent HEGs having complex optical morphologies (no upper limits on 
the optical core luminosities could be estimated). 
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2 4 6 8 10 




Fig. 6. Same as Fig.0 with over plotted the 'lines' Li on /vL^'\ oha as function of Lo / L-K'\ ohs: obtained from eq. Eland 
01 for several values of / and Ay. torus- The dot-dashed lines correspond to fixed value of / and running Ay, torus, and 
vice- versa for the dotted lines (as labelled) . 
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Fig.D.l. HST/WFPC2 F702W images (left-hand panels) and ARNICA/NICS images (right-hand panels) of the 
HEG-FRII galaxies of the sample. 
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Fig.D.2. (cont.) 
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